Chapter 3 



Theory of Quadratic TFDs 



The quadratic time-frequency distributions (TFDs) introduced in the last chapter 
represent the majority of the methods used in practical applications that deal with 
non-stationary signals. In this chapter, which completes the introductory tutorial, 
we show that these particular quadratic TFDs belong to a general class of TFDs 
whose design follows a common procedure, and whose properties are governed by 
common laws. This quadratic class may be considered as the class of smoothed 
Wigner-Ville distributions (WVDs), where the “smoothing” is described in the ( t , /) 
domain by convolution with a “time- frequency kernel” function y(t, /), and in other 
domains by multiplication and/or convolution with various transforms of 7 (t,f). 
The generalized approach allows the definition of new TFDs that are better adapted 
to particular signal types, using a simple and systematic procedure as opposed to 
the ad hoc methods of Chapter 2. 

The first section (3.1) extends Section 2.1 by enumerating in detail the key 
properties and limitations of the WVD. Thus it motivates the introduction of general 
quadratic TFDs (Section 3.2) and prepares for the discussion of their properties 
(Section 3.3). In Section 3.2, using Fourier transforms from lag to frequency and 
from time to Doppler (frequency shift), the quadratic TFDs and their kernels are 
formulated in four different but related two-dimensional domains. One of these, 
namely the Doppler-lag ( v , r) domain, leads to the definition of the “ambiguity 
function” and allows the smoothing of the WVD to be understood as a filtering 
operation. In Section 3.3, the list of properties of the WVD is supplemented by 
mentioning some desirable TFD properties not shared by the WVD. The various 
TFD properties are then expressed in terms of constraints on the kernel, so that 
TFD design is reduced to kernel design. Three tables are provided showing the 
kernel properties equivalent to various TFD properties, the kernels of numerous 
popular TFDs in the various two-dimensional domains, and the properties of those 
same TFDs. 



0 Author: Boualem Boashash, Signal Processing Research Centre, Queensland University 
of Technology, Brisbane, Australia. Reviewers: K. Abed-Meraim, A. Beghdadi, M. Mesbah, 
G. Putland and V. Sucic. 
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3.1 The WVD 

3.1.1 Properties of the WVD 

We have seen that the WVD has the simplest time-lag kernel (see Eq. (2.7.7)), and 
that each of the other TFDs can be written as a filtered WVD using a specific time- 
lag kernel filter (see Eq. (2.7.37)). In this sense the WVD is the basic or prototype 
TFD and the other TFDs are variations thereon. Moreover, we see in Fig. 2.7.1 that 
the WVD gives the sharpest indication of the IF law of a linear FM signal. These 
are some of the reasons why the WVD is the most widely studied TFD and deserves 
a further detailed description of its properties, as listed below. 

• Realness (RE): W z (t. /) is real for all z, t and /. 

• Time-shift invariance (also called time covariance): A time shift in the 
signal causes the same time shift in the WVD; that is, if 

z r (t) = z(t - t 0 ) , (3.1.1) 

then 

W Zr (t,f) = W z (t-t 0 ,f). (3.1.2) 

• Frequency-shift invariance (also called frequency covariance): A fre- 
quency shift in the signal causes the same frequency shift in the WVD; that 
is, if 

z r (t) = z(t) e j27r ^ ot , (3.1.3) 

then 

W Xr (t,f) = W z (t,f-f 0 ). (3.1.4) 

• Time marginal (TM): Integration of the WVD over frequency gives the 
instantaneous power: 

/ OO 

W z (t,f)df = \z(t)\ 2 . (3.1.5) 

-oo 

• Frequency marginal (FM): Integration of the WVD over time gives the 
energy spectrum: 

/ OO 

W z (t,f)dt=\Z(f)\ 2 . (3.1.6) 

-oo 



• Global energy: Integration of the WVD over the entire (t, /) plane yields 
the signal energy E z : 



/»00 nOO 

/ / W z (t,f)dtdf = E z . 



— oo J — OO 



(3.1.7) 
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• Instantaneous frequency (IF): For an analytic signal, the first moment 
[i.e. the mean] of the WVD w.r.t. frequency is the IF: 



JZof W ^f)df 



1 d r , 

s:s [ “ g2(i)1 



(318) 



• Time delay (TD): The first moment of the WVD w.r.t. time is the TD: 



Qw*(*./)<** 

SZo W ^f)dt 



-TW l “ gZIJ)h 



(3.1.9) 



where Z(f) is the FT of z(t). 

• Time support (TS): The time support of W z (t, /) is limited by the duration 
of z{t)\ that is, if z(t) = 0 for t < t\ and for t > t 2 , then W z {t, /) = 0 for 
t <t\ and for t > t 2 . 



• Frequency support (FS): The frequency support of W z (t,f) is limited by 
the bandwidth of z{t)\ that is, if tF{z(t)} = Z(f) = 0 for f < f\ and for 
/ > h , then W z (t, /) = 0 for / < /i and for / > / 2 . 

• Convolution invariance: The WVD of the time-convolution of two signals 
is the time- convolution of the WVDs of the two signals; that is, if 



Z 3 (t) = zx(t) * z 2 (t) , (3.1.10) 

then 

W Za (t,f) = W Zl (t,f)*W Xa (t,f). (3.1.11) 



• Modulation invariance: The WVD of the frequency- convolution of two 
signals is the frequency- convolution of the WVDs of the two signals; that is, 
if 

Z 3 (t) = zi(t)z 2 (t) , (3.1.12) 



then 



W Z3 (t,f) = W Zl (t,f)*W Z2 (t,f). 



(3.1.13) 






Invertibility: If W z (t. f) is the WVD of the signal z(t), it may be shown [1, 
pp. 223-4] that 



W z (t/2J)e^ t df = z(t)z*(0). 



(3.1.14) 



Putting t = 0 in this result yields the magnitude, but not the phase, of z(0). 
Hence a signal z(t) may be recovered from its WVD up to a phase factor. 
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• Inner-product invariance: The WVD is a unitary transformation; that 
is, it preserves inner products: 




W z i (t, f ) W Z2 (t, f) dt df 



Zi(t) z^{t) dt 



(3.1.15) 



The above properties are not independent; for example, it is a trivial exercise to 
show that either of the marginals (TM or FM) implies the global energy condition. 

A more comprehensive list of the WVD properties is given in [1] . These either 
follow directly from the definition, or are proven in [2]. 

A notable omission from the above list is positivity; the WVD can assume 
negative values, and indeed does so for almost every signal. 

Applicability: Many properties of the WVD are desirable in applications; for 
example, realness is consistent with the notion of energy density in as much as 
energy is real, while convolution invariance and modulation invariance make the 
WVD partly compatible with linear filtering theory. Nevertheless, the WVD was 
not applied to real-life problems until the late 1970s, when it was implemented for 
the purposing of processing linear FM signals used in seismic prospecting [3] . 



3.1.2 Limitations of the WVD 

Despite its many desirable properties, the WVD has some drawbacks. It may assume 
large negative values. Further, because it is bilinear in the signal rather than linear, 
it suffers from spurious features called artifacts or cross-terms, which appear 
midway between true signal components in the case of multicomponent signals as 
well as non-linear mono- and multicomponent FM signals [see Article 4.2]. 



3. 1.2.1 Nonlinear Monocomponent FM Signals and Inner Artifacts 

Consider a monocomponent signal z(t) = a(t ) e^u) . In the case of a linear FM 
signal, the WVD gives an accurate representation of the IF law (Fig. 1.1.3) because 
the CFD approximation to is exact, so that the signal kernel is a dechirped 
function of r [see Eqs. (2.1.19) to (2.1.25)]. In the case of a nonlinear FM signal, the 
CFD approximation is not exact and the signal kernel (i.e. the IAF) is not dechirped. 
This results in the formation of inner artifacts, which arise from “within” a single 
component. 

An example of a nonlinear FM signal is the hyperbolic FM signal with rect- 
angular amplitude a(t) and the phase given by 



<j>(t) = In |1 + at|. 



The IF is 



fi(t) = 



y(t) 



fo 



2n 1 + at 

For the time interval 0 < t < T , the starting frequency is fo 
frequency is f max = fi(T) = f 0 /(l+aT) . 



(3.1.16) 

(3.1.17) 
fi( 0). The finishing 
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Fs=1Hz N=65 Fs=1Hz N=65 

Time-res=1 Tlme-res=1 




(a) (b) 



Fig. 3.1.1: WVD of (a) hyperbolic FM signal with starting frequency 0.1 and finishing frequency 
0.4 Hz; (b) linear FM signal with starting frequency 0.4 and finishing frequency 0.1 Hz. Both plots are 
for a duration of 65 samples (sampling rate 1Hz). 



Fig. 3.1.1(a) shows the WVD of a hyperbolic FM signal with duration 65 samples, 
fo = 0.1 and f ma x = 0.4. While the crest of the WVD seems to be a reasonable 
approximation to the IF law, the energy concentration is poorer. The many spurious 
ridges are the inner artifacts. They alternate in sign as we move normal to the IF 
law in the (f, /) plane; this is a characteristic feature of inner artifacts. 

For comparison, Fig. 3.1.1(b) shows the WVD of the linear FM signal with the 
same duration and the same frequency limits (with falling frequency). Note the 
superior energy concentration for the linear FM case, and the more attenuated 
artifacts distributed on both sides of the main ridge. 

Artifacts caused by nonlinear FM laws can be reduced by windowing the IAF in 
r before taking the FT, leading to the windowed WVD. This procedure, however, 
causes a loss of frequency resolution [4]. 

3. 1.2. 2 Multicomponent Signals and Outer Artifacts 

If z(t) is a multicomponent signal, the algebraic expansion of K z (t,r) contains 
cross-product terms which, when Fourier-transformed, give rise to spurious features 
in the WVD. These are the outer artifacts or cross-terms. To explain these, 
consider the signal 

z(t) = Zi(t) + 22 (f) (3.1.18) 

where z(t), z\(t) and 22 (f) are analytic. Expanding the IAF, we obtain 

K z (t,r) = K Zl (t,r) + K Z2 (t,r) + K ZlZ2 (t,r) + K Z2Zl (t,r) (3.1.19) 

where K ZlZ2 (t,r) and K Z2Zl (t,r) are the “signal cross-kernels” or instantaneous 
cross-correlation functions (e.g. K ZlZ2 (t , r) = Z\(t + ^) 2 jj(f — |) ). Taking FTs of 
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Fs=1 Hz N=65 
Time-res=1 
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Fig. 3.1.2: The WVD of the sum of two linear FM signals with frequency ranges 0.1— 0.2 and 0.3— 0.4 
Hz, unit amplitudes, and a duration of 65 samples (sampling rate 1Hz). 



Eq. (3.1.19) w.r.t. r [using Eq. (2.1.18)], we find 

w z (t, f ) = W Z1 ( t , /) + W Z2 ( t , /) + 2Re{W ZlZ2 ( t , /)} (3.1.20) 

where W Zl (t, /) and W Z2 (t, f) are the WVDs of Zi(t) and Z 2 (t), and W ZlZ2 (t, f) is 
the cross-Wigner-Ville distribution (XWVD) of Zi(t) and defined by 

/ OO 

Zl (t+l)z* 2 (t-i)e-^ T dT. (3.1.21) 

-OO 

Thus the WVD of the sum of two signals is not just the sum of the signals’ WVDs, 
but also of their XWVDs. If Z\{t) and Z 2 (t) are monocomponent signals, W Zl (t. f ) 
and W Z2 (t,f) are the auto-terms, while 2Re{IV ZlZ2 (f, /)} is a cross-term. 

Fig. 3.1.2 shows the WVD of the sum of two parallel linear FM signals. There 
seem to be three components rather than two; the “extra” component at the mean 
frequency of the expected components has a large oscillating amplitude, and occurs 
in a region of the (f, /) plane where we expect no energy at all. Fig. 3.1.3(a) shows 
the WVD of the sum of two FM signals crossing over in frequency. A large number 
of undulations appear in addition to the two main ridges. In both examples, the 
cross-terms alternate in sign as we move parallel to the expected features in the 
(t, f) plane; this is a characteristic feature of cross-terms. 

3. 1.2. 3 Suppression of Cross-Terms 

Cross-terms can make the WVD difficult to interpret, especially if the components 
are numerous or close to each other, and the more so in the presence of noise; cross- 
terms between signal components and noise exaggerate the effects of noise and 
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Fig. 3.1.3: Suppression of cross-terms in the sum of a rising hyperbolic FM signal (frequency range 
0.1 to 0.4 Hz) and a falling linear FM signal (frequency range 0.4 to 0.1Hz), with unit amplitudes and 
a duration of 65 samples (sampling rate 1Hz): (a) WVD; (b) spectrogram with 21-point rectangular 
window; (c) masked WVD, being the product of (a) and (b). 



cause rapid degradation of performance as the SNR decreases. For such reasons, 
cross-terms are often regarded as the fundamental limitation on the applicability 
of quadratic time-frequency methods, and the desire to suppress them has led to 
several approaches such as: 

1. If we multiply the WVD by the spectrogram, we obtain the so-called masked 
WVD [5] , which combines the cross-term suppression of the spectrogram with 
the high resolution of the WVD. To illustrate the effect, Fig. 3.1.3 shows the 
WVD, spectrogram and masked WVD of the sum of two crossed FM signals; 
observe that the masked WVD is “cleaner” than the WVD and has better 
resolution than the spectrogram. 

2. Eq. (3.1.21) indicates that the XWVD of Z\(t) and 22 (f) is linear in 22 (f). If 
21 (f) is a reference signal and 22 (f) is the signal under analysis, there will be 
no cross-terms between components of 22 (f)- This observation inspired efforts 
to displace the WVD with the XWVD in relevant areas of application. A 
reference signal is readily available in optimal detection [6], sonar and radar [7], 
and seismic exploration [8]. In applications where reference signals are not 
available, a filtered version of the observed signal is used as a reference signal. 
The filtering procedure may use the IF as a critical feature of the signal; for 
example, the data-dependent TFDs defined in [9] use, as a reference signal, 
the signal component that maximizes the energy concentration in the TFD. 

3. All quadratic TFDs in Chapter 2 can be identified by a time-lag kernel, includ- 
ing the spectrogram which suppresses cross-terms. By analyzing the properties 
of time-lag kernels, we can define and design TFDs that attenuate cross-terms. 
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A quadratic TFD in which the cross-terms are attenuated relative to the auto- 
terms is often called a reduced-interference distribution (RID). 

This chapter concentrates on the RID approach for defining quadratic TFDs with 
the property of cross-terms suppression. 



3.2 Formulations of Quadratic TFDs 

3.2.1 Time-Lag Formulations and Other Domain Definitions 

Given an analytic signal z(t), the instantaneous autocorrelation function (IAF) was 
defined in Eq. (2.1.8) as 



K z (t,r) = z(t + l)z*(t - §) (3.2.1) 

The IAF K z (t,r) is a function of 2 time variables: the actual time t and the time 
lag t. By taking the dual domain of t and t in frequency, we obtain the frequency 
variables v and /. This allows for four domains of representation: 

(t,r) 

(tJ) 

(g/) 

(gt). 

Starting with the IAF, we define the Wigner-Ville distribution (WVD) by taking 
the FT (r — > /): 

W x (t,f)= T{K x {t,r)}. (3.2.2) 

T — ►/ 

The FT ( t — > v) of the WVD defines the spectral autocorrelation function (SAF) 
[Eq. (2.1.26)] as 

k x (v,f) = £ u {W x (t,f)}. (3.2.3) 

The FT (t — > v) of the IAF K z (t,r) equals the IFT (r <— /) of the SAF k z {v,f ), 
and defines the symmetrical ambiguity function (AF) as: 

A x {v,t) = T{K z {t,T)} = {Mi/,/)} (3-2-4) 

where . .} represents the forward (direct) Fourier transform, and . .} its 

inverse. 1 

1 Many authors define the symmetrical ambiguity function using an inverse transform from t 
to is, with the result that the arrows in the “diamond diagrams” [e.g. Eq. (3.2.5)] point upwards 
instead of to the right (this affects only the lower- left and upper-right arrows). Here we choose 
our definitions so that transforms from a time variable ( t or r) to a frequency variable (is or /) 
are always forward. 
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If we represent Fourier transformations by arrows labeled with the participating 
variables, Eqs. (3.2.2) to (3.2.4) may be combined into the single graphical equation 

W z (tJ ) 




K z (t,r) 




A z {v,t) 



k z {v,f). 




(3.2.5) 



This single graphical representation of several equations is very useful in that it links 
several known methods by simple FT operations. The knowledge of FT properties 
then allows to use characteristics of a method in one domain and transfer them 
to another domain. We can so relate radar methods in the Doppler-lag domain to 
(t, /) methods and to cyclostationary methods. 

Now let us represent similarly the TFD kernel in these four domains by taking 
various FTs of the time-lag kernel G(t, r): 



7 (t,f) 




(3.2.6) 



We refer to g{y,r) as the Doppler-lag kernel, to Q(y,f) as the Doppler- 
frequency kernel, and to j(t, /) as the time-frequency kernel. 



3.2.2 Time-Frequency Formulation 

We have defined the smoothed IAF as 

R z (t,T)=G(t,T)*K z (t,r) (3.2.7) 

where G(t,r) is the time-lag kernel. Following Eq. (2.7.1), we define a class of 
quadratic TFDs (which are smoothed WVDs) as 

p z (t,f)= f{R z (t,r)}. (3.2.8) 

T— ►/ 

Substituting Eq. (3.2.7) into Eq. (3.2.8), writing out the convolution and the trans- 
form, and substituting from Eq. (3.2.1), we obtain 

/ OO /»oo 

/ G(t—u,T)z(u+^)z*(u — : ^)e~ :i2 ' KfT dudT (3.2.9) 

-ooJ—oo 




68 



Chapter 3: Theory of Quadratic TFDs 



which defines the general class of quadratic TFDs in terms of the signal and the 
time-lag kernel. The class of TFDs of this form is called the quadratic class . 2 
Eq. (3.2.8) is included in the graphical equation 

Pz(t,f) 




which assigns symbols to the various FTs of the smoothed IAF [11, p.436]. By 
analogy with Eqs. (3.2.3) and (3.2.4), we call r z (v,f) the “generalized SAF” and 
A z (v,t) the “generalized ambiguity function” (GAF) 3 or more precisely filtered 
ambiguity function; note the distinction between the normal A in Eq. (3.2.5) and 
the calligraphic A in Eq. (3.2.10). 

Using Eq. (3.2.7) and the lower two arrows in Eq. (3.2.10), together with the 
convolution properties of the FT, we obtain in sequence 

A z (u,t) = g(v,r) A z (v,t) (3.2.11) 

‘r z {v,f) = G{v,f)*k z {v,f) (3.2.12) 

Then, using Eq. (3.2.7) and the upper left arrow in Eq. (3.2.10), together with the 
convolution properties, we obtain 

Pz(t, f) = l{t, /)** W z (t,f) (3.2.13) 

where the double asterisk denotes double convolution in t and / (cf. [11], 
pp. 437,475). 

Eq. (3.2.13) defines the quadratic class of TFDs in terms of the WVD and the 
time-frequency kernel. For this reason we regard the WVD as the “basic” or “pro- 
totype” quadratic TFD, and all other quadratic TFDs as filtered versions thereof. 



2 The quadratic class as defined here satisfies the time-shift- and frequency-shift-invariance 
properties. Most current authors seem to equate “Cohen’s class” with the quadratic class. In 
addition, some authors apply the term “quadratic time-frequency representation” or “QTFR” to 
all time-frequency representations that are quadratic in the signal (whether they have the form of 
Eq. (3.2.9) or not) so that “Cohen’s class” becomes a subset of the “QTFRs”. However, the class 
originally defined by Cohen in a quantum-mechanical context [10] differs from the quadratic class 
[Eq. (3.2.9)] in that the marginal conditions must be satisfied, while the kernel may depend on the 
signal (in which case the TFD is not quadratic in the signal). We shall see (e.g. in Table 3.3.3) 
that many useful quadratic TFDs do not satisfy the marginals, so that they are not members of 
“Cohen’s class” as originally defined. 

3 The term generalized ambiguity function is also used with a different meaning in connec- 
tion with polynomial TFDs; see Article 5.5. 
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Eq. (3.2.13) also suggests a method of suppressing cross-terms in the WVD. Be- 
cause inner artifacts alternate as we move in the frequency direction, they can be 
attenuated by choosing a 7 (t, /) with a sufficient spread in the frequency direction. 
Similarly, because cross-terms (outer artifacts) alternate as we move in the time 
direction, they can be attenuated by choosing a j(t, f ) with a sufficient spread in 
the time direction. 

Eqs. (3.2.10) and (3.2.11) show how one quadratic TFD can be transformed to 
another by a linear transformation. 

3.2.3 Doppler-Lag Formulation and TFD Design 

Substituting for K z (t, r) in Eq. (3.2.4) and writing out the transform, we obtain 

/ OO 

z{t+^)z*(t-^)e~ i2KVt dt (3.2.14) 

-OO 

/ OO 

z(t + §) [z{t - § ) e j2nvt ] * dt. (3.2.15) 

-OO 

The expression in square brackets can be obtained by delaying z(t + ■§) in time by 
r and shifting it in frequency by v, indicating that A z r) is the correlation of the 
signal with a time-delayed and frequency-shifted version of itself. This correlation 
is well known in radar theory as the Sussman ambiguity function [12]; the name 
“ambiguity” arises from the equivalence between time-shifting and frequency-shifting 
for linear FM signals, which are frequently used in radar. Hence the Doppler-lag 
(17 r) domain is also called the ambiguity domain. 

Eq. (3.2.10) indicates that the quadratic TFD is a 2D FT (half inverse, half 
forward) of the filtered ambiguity function. Writing out the transforms gives 

/ OO /»oo 

/ g(v,T)A z (v,T)e j2n{l ' t -- fT) dvdT. (3.2.16) 

-00J—00 

Then writing the dummy u for t in Eq. (3.2.14) and substituting the result into 
Eq. (3.2.16) gives 

p z {t,f) = jjj g{iy,T)z{u+^)z*{u-^)e j2 ^ vt - vu - fT) dududT (3.2.17) 

where the integrals are from —00 to 00 ; this defines the quadratic TFD in terms of 
the Doppler-lag kernel 

From the lower left side of Eq. (3.2.6), the relationship between the time-lag 
kernel and the Doppler-lag kernel is 

g(v,T) = £ v {G(t,T)}. (3.2.18) 

Using this equation, the time-lag kernels determined in Section 2.7 may be converted 
to Doppler-lag form. For example, using Eq. (2.7.7), we find that the Doppler-lag 
kernel for the WVD is g{v,r) = T {<$(£)} = 1 . 
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Eq. (3.2.13) and its 2D FT Eq. (3.2.16) express that a quadratic TFD can be 
designed using basic filter design principles. As in ID filter design, the filter specifi- 
cations are best given in the domain where the filtering operation is expressed as a 
multiplication as in Eq. (3.2.11), rather than in the dual domain with a convolution 
as in Eq. (3.2.13). 

If the Doppler-lag kernel g(u, r) has the separable form G\(v) ^(t), multiplica- 
tion by this kernel may include the combined effect of time windowing and frequency 
windowing (see the discussion of Eq. (2.1.40)). 

We can use the same process to define a filter which attenuates cross-terms in 
quadratic TFDs. The cross-terms in the (i, /) domain tend to be highly oscillatory, 
so that the corresponding terms in the dual (;/, r) domain tend to be far from 
the origin (high- pass). The auto-terms in the (t, f) domain tend to be smooth 
and well delineated, so that the corresponding terms in the dual {v, r) domain are 
concentrated about the origin or “pass through” the origin (low-pass) . This behavior 
is well known in the field of radar [13, 14]. Hence the cross-terms in the ambiguity 
domain can be “filtered out” by selecting a kernel filter 4 g(z/, r) that deemplrasizes 
information far from the origin in the Doppler-lag domain. 

Various authors [10, 15, 16] have shown that other desirable TFD properties are 
equivalent to constraints on the kernel filter, and we shall see that most of these 
constraints are conveniently expressed in the Doppler-lag domain. For all the above 
reasons, the design of the TFD kernel filter is usually performed in the Doppler-lag 
domain. Often, the resulting kernel is then described in the time-lag domain for 
ease of TFD implementation, as given by Eq. (3.2.9). 

3.2.4 Doppler-Frequency Formulation 

From Eqs. (3.2.10) and (3.2.12) we have 

Pz(t, f) = T~l {< 3{v , /) * k x (u, /)} (3.2.19) 

or, writing out the transform and convolution, 

/ OO /»oo 

/ Q{i>,f~il)k z {u,'n)e :,2 ' KVt dr]du (3.2.20) 

-OO*' — OO 

Writing 77 for / in Eq. (2.1.31) and substituting the result into Eq. (3.2.20), we obtain 

/ OO /*00 

/ g(v,f-v)Z( V +%)Z*(r 1 -%)ei 2 ™ t d V dv. (3.2.21) 

-00*/ — OO 

This defines quadratic TFDs in terms of the Doppler- frequency kernel f) and 
the signal spectrum. The interest of this formulation is that TFDs of narrow-band 
signals expressed by their spectra may be more efficiently computed in this form. 

4 The term “kernel filter” is used to reinforce the idea that designing quadratic TFDs essentially 
reduces to filter design with specific constraints. 
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3.2.5 Examples of Simple TFD Formulations 
3. 2. 5.1 Separable Kernels 

A simple way to design kernel filters for quadratic TFDs is 


to consider the special 


case of a separable kernel (see Article 5.7): 

g(”,T) = Gi(v) g 2 (r). 


(3.2.22) 


If we let 

Gi{v) = £{gi{t)} 


(3.2.23) 


and 

G 2 (/)= F {g 2 (r)}, 

T—>f 


(3.2.24) 


then Eq. (3.2.11) becomes 

A z (v,t) = Gi (v) g 2 (t) A, (v,t). 


(3.2.25) 


Then, using Eqs. (3. 2. 10), (3.2.23), and(3.2.24), we obtain 
Pz(t,f) = gi(t) *W z (t, f) *G 2 (f). 


(3.2.26) 



Eq. (3.2.25) shows that the design of the kernel filter G\{v) g 2 (r) is greatly sim- 
plified as the 2D filtering operation is replaced by two successive ID filtering op- 
erations. Equivalently, in Eq. (3.2.26), the two convolutions can be evaluated in 
either order, indicating that the Doppler-dependent and lag-dependent factors in 
the separable kernel g(u 1 r) lead to separate convolutions in time and frequency, 
respectively. 

A separable kernel is not to be confused with a product kernel, which is a 
function of the product vt , e.g. 

9{v,t)= gs{vT). (3.2.27) 

3. 2. 5. 2 Doppler-Independent Kernels 

A Doppler-independent (DI) kernel is a special case of a separable kernel ob- 
tained by putting 

G\{v) = 1 (3.2.28) 

in Eqs. (3.2.22) and (3.2.23), which then become 

9{v, T ) = 52 (t) (3.2.29) 

gi{t) = S{t). (3.2.30) 

Making these substitutions in Eqs. (3.2.25) and (3.2.26), we obtain 

A z {v,t) = g 2 {r) A z {v,t) (3.2.31) 

p z (t,f) = G 2 (f)*W z (t,f). (3.2.32) 

The last result shows that a DI kernel is obtained by applying only ID filtering and 
causes smearing of the WVD in the frequency direction only. 
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3. 2. 5. 3 Lag-Independent Kernels 

A lag-independent (LI) kernel is another special case of the separable kernel, 
obtained by putting 

92(t) = 1 (3.2.33) 

in Eqs. (3.2.22) and (3.2.24), which then become 

g(v,r) = Gi(u) (3.2.34) 

G 2 (f) = 6(f). (3.2.35) 

Making these substitutions in Eqs. (3.2.25) and (3.2.26), we obtain 

A z (u, t ) = Gi(v) A z (y, t) (3.2.36) 

Pz(t,f) = 9i(t)*W z (t,f). (3.2.37) 

The last result shows that an LI kernel is obtained by applying only ID filtering 
and causes smearing of the WVD in the time direction only. 

The kernel of the WVD is g(v,r) = 1, which is both Doppler-independent and 
lag-independent. The windowed WVD kernel, however, is Doppler-independent. 

Article 5.7 provides an in-depth treatment of such separable kernels, including 
examples of quadratic TFDs obtained by this simple kernel filter design procedure. 

3.3 Properties of Quadratic TFDs 

3.3.1 Desirable Properties 

Some quadratic TFDs verify desirable properties that are not shared by the WVD, 
and vice versa. Later we shall relate the properties of a TFD to the constraints 
of its kernel, and tabulate properties for selected TFDs. But first we discuss some 
properties that have been promoted [17-19] as fundamental for a wide range of 
applications. 

1. Concentration of local energy: The energy in a certain region R in the 
(t,f) plane, denoted by E ZR , is given by the integral of the TFD over the 
region R; for example, if R is the region within a time interval At and a 
frequency band A/, the energy within R is 

E zr = f [ p z (t,f)dfdt. (3.3.1) 

JAt JAf 

2. IF/TD visualization: The TFD of a monocomponent signal directly depicts 
the instantaneous frequency /,(t) and time delay ra(/) as a range of peaks 
along the curve representing the FM law in the (t, /) plane. That is, if z(t) is 
a monocomponent signal and p z (t, f) is its TFD, then 

max (i, /) = fi(t). 



(3.3.2) 
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3. Reduced interference (RI): The TFD attenuates or suppresses inner and 
outer artifacts (cross-terms) relative to the signal components (auto-terms) . 

Another condition considered by some authors is positivity, or more precisely 

non-negativity (NN), defined as: 

Pz(t,f)> 0 V z,t,f. (3.3.3) 

Among quadratic TFDs, only sums of spectrograms possess NN. This means that 
the Doppler-lag kernel of a “positive” TFD is a sum of ambiguity functions, which 
makes the NN property incompatible with both the IF property (see Table 3.3.1) 
and IF visualization [20] . Hence NN is usually considered non-essential because its 
cost is excessive. 

A subclass V of quadratic TFDs comprises those TFDs which satisfy the re- 
alness, time marginal, frequency marginal, instantaneous frequency, time support 
and frequency support properties. Researchers have shown much interest in TFDs 
of this class, and especially RIDs of this class. To some extent the design of RIDs 
is the art of “improving” on the spectrogram by 

• sacrificing NN, 

• improving resolution, and 

• retaining sufficient reduced-interference ability for the application. 

The result is usually a compromise between the spectrogram and the WVD, involv- 
ing a time-frequency kernel filter less extensive than that of the spectrogram. The 
compromise may involve sacrificing one or both of the marginal properties (time or 
frequency). Although the marginals are critical in the field in quantum physics [10], 
they seem to be less important in signal processing. As evidence of this we may cite 
the following: 

• The spectrogram does not satisfy the marginals and yet has always been re- 
garded as a very useful TFD. 

• Attempts to improve on the spectrogram are most often motivated by the 
need for higher resolution (see e.g. [21]) rather than any desire to satisfy the 
marginals. 

• If the marginal conditions do not hold, it is still possible that the integral forms 
thereof [Eqs. (2.6.6) and (2.6.7)] are approximately true over sufficiently broad 
time intervals or frequency bands. 

For these reasons it is suggested that V is not necessarily appropriate for signal- 
processing applications, and an alternative class V , which comprises reduced- 
interference distributions satisfying the realness, global and local energy, IF vi- 
sualization and components resolution properties, would be more relevant. Note 
that the properties of class V are the same as the properties listed at the end of 
Section 1.1.5. 
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3.3.2 TFD Properties & Equivalent Kernel Constraints 

Using some properties listed in [22], the relevant kernel constraints were adapted 
for separable, DI and LI kernels. The results are collected in Table 3.3.1. 

The proofs of these kernel constraints use Eq. (3.2.13), which states that the 
quadratic TFD is the WVD convolved in t and / with the time-frequency kernel 

7 

Pz(t, f) = 7 (t, f ) ** W z (t, /). 

Using this relation, we can find siLfficient conditions under which certain properties 
of the WVD carry over to p z (t, /). For example: 

• Realness holds if 7 (t,f) is real. 

• Because 2D convolution is shift-invariant, time- and frequency-shift invariance 
hold for any fixed 7 (t, /). 

• Time support and time extent hold for a DI kernel, which does not redistribute 
the WVD in time, whereas frequency support and frequency extent hold for 
an LI kernel, which does not redistribute the WVD in frequency. 

• The IF moment property holds for a DI kernel if G* 2 (/) has a first moment 
(mean frequency) of zero, so that convolution w.r.t. / does not change the 
first moment of the WVD w.r.t. /. Similarly, the GD property holds for an 
LI kernel if gi{t) has a first moment (mean time) of zero, so that convolution 
w.r.t. t does not change the first moment of the WVD w.r.t. t. 

Some further proofs of kernel constraints are given by Cohen [23] . 

3.3.3 Examples of TFDs with Specific Properties 

For the TFDs defined so far, Table 3.3.2 lists their kernels in various domains, and 
Table 3.3.3 shows their properties. 

Note: A time-lag kernel satisfying the time-support constraint is described as a 
butterfly function [11] or cone-shaped kernel [24]; that is, the nonzero values 
of the kernel are confined to the interior of a two-dimensional cone in the (f , r) plane 
(see the entry in the “General” column of Table 3.3.1). The Born- Jordan and ZAM 
distributions have kernels of this type. 

Inspection of Table 3.3.1 confirms that the WVD satisfies all of the listed prop- 
erties except reduced interference (ri). The reduced-interference capabilities of sep- 
arable kernels warrant special attention. Because the inner artifacts alternate as we 
move normal to the components, they also alternate in the frequency direction and 
can be suppressed by convolution with a sufficiently long G 2 {f) [see Eq. (3.2.32)], 
which corresponds to a sufficiently short <72 (t). This is possible for a DI kernel but 
not an LI kernel. Because the cross-terms (outer artifacts) alternate as we move 
parallel to the components, they also alternate in the time direction and can be 
suppressed by convolution with a sufficiently long gi(t) [see Eq. (3.2.37)], which cor- 
responds to a sufficiently short G\(y). This is possible for an LI kernel but not a 
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Table 3.3.1: TFD properties and associated constraints on the Doppler-lag kernel for general, sep- 
arable, Doppler-independent and lag-independent kernels. The asterisk in “WVD only*” means that 
the WVD may be multiplied by a scale factor. Abbreviations for properties: NN = non-negativity; 
RE = realness; TI = time-shift invariance (or time covariance); FI = frequency-shift invariance (or fre- 
quency covariance); TM = time marginal; FM = frequency marginal; IF = instantaneous frequency; 
TD = time delay; TS = time support; FS = frequency support; RI=RID capability. 





KERNEL CONSTRAINTS 


Property 


General 


Separable 

g(v,r) = G\(v) 32 (r) 


DI 

S(J'T) = 52 (r ) 


LI 

g{v,r) = Gi(u) 


NN: 


G(t, t ) is a 
sum of IAFs. 


5i(t)fl2(r) is a 
sum of IAFs. 


Never 


Never 


RE: 


fl(G-r) 

= g*(-v, -r). 


Gi{v) 52 (r) 

= G\(-v) g^-T). 


G2(/) is real. 


51 (t) is real. 


TI: 


g{v,r) does not 
depend on t. 


G\{v) </2(r) does 
not depend on t. 


<72 (r) does not 
depend on t. 


Gi(i/) does not 
depend on t. 


FI: 


g(v , r) does not 
depend on /. 


Gi(y) g2(r) does 
not depend on /. 


52 (r) does not 
depend on /. 


Gi(i/) does not 
depend on /. 


TM: 


g{v, 0) = 1 Mv 


Gi{u) 32(0) = 1 V12 


52 (0) = 1 


WVD only 


FM: 


5(0, r) = 1 Vr 


Gi( 0 ) p 2 (r) = 1 Vr 


WVD only 


Gi( 0 ) = 1 


IF: 


g(v, 0) = const. 

!fUo=0 VzA 


G i(i>) 52(0) = const. 
52 ( 0 ) = 0. 


52(0) = 0 


WVD only* 


TD: 


5(0, r) = const. 
fUo =0 Vr. 


Gi( 0 ) 52 (r) = const. 
G'i ( 0 ) = 0 . 


WVD only* 


O 

II 

O 


TS: 


G(i,r) = 0 
if |r| < 2|t|. 


DI only 


Always 


WVD only* 


FS: 


S(f,v) = 0 

if \u\< 2 \f\. 


LI only 


WVD only* 


Always 


RI: 


Unrestricted 


Unrestricted 


Inner x-terms 


Outer x-terms 



DI kernel. A general separable kernel with sufficiently short G i(z') and (j'lij) can 
therefore attenuate both kinds of artifacts; a TFD with such a kernel is shown in 
Fig. 2 . 7.1 part ( 1 ) for a linear FM signal. 

Properties of separable kernels and their special cases, such as the B-distribution 
and Modified B-distribution, are discussed further and illustrated on examples in 
Article 5 . 7 . Design of RIDs is discussed more generally in Article 5 . 2 . 
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Table 3.3.2: Kernels of selected TFDs in the time-lag, Doppler-lag, and (where possible) Doppler- 
frequency domains. The window w(t) is assumed to be real and even. Its FT and AF are W (/) and 
A w (iy,r), respectively. The prefix “w-" means “windowed”. 



Distribution 


KERNEL 


G(t,r) 


s(D' r) 


Qiy,f) 


Wigner-Ville 


m 


1 


5(f) 


Levin 


i [*(t+§) + *(*-§)] 


cos (tvi/t) 


hW+?) + «(f- 1)] 


Born-Jordan 


1357T rect 537 


sine (2 out) 


i rert J— 

\2cv\ lav 


Modified B 


cosh - t 


|r(/3+j7n/)| 2 


\T{P+j-nv)\ 2 


f 3°^ cosh- 2 P £ d£ 


r*W) 


T^{p) °\J > 


tu-WVD 


S(t) w(t) 


w(t ) 


W(f ) 


w-Levin 


^l[S(t+^ + 5(t-^\ 


w(t) cos (tyvt) 


h[w(f+%) + w{f-%)] 


ZAM 


W ( T ) reCt 2 7JZ 


wM^f-rsinc^ 
k ' 2|t| a 




Rihaczek 


S(t- 1) 


„ — ITT 1ST 

e J 


<5(/+|) 


w-Rihaczek 


w{T)S{t- |) 


w(r) e~^ VT 


W(f+ |) 


Page 


S{t- III) 


e ~j-nv\T\ 


+ ju/[2n(f 2 -u 2 /4:)} 


Choi-Williams 


y/Ta — 7r 2 crt 2 /r 2 

M e 


e T /<7 


VTa -tt 2 a f 2 /u 2 

W 1 


B 


\t\P cosh _2/3 1 


hi' 3 \r(p+jirv)\ 2 

21-2/3 T(2/3) 




Spectrogram 


w(t+ §)w(t- §) 


A w (v,t) 


W(f+$)W{f-%) 



3.4 Summary, Discussion and Conclusions 

This chapter ends the tutorial introduction to TFSAP constituted by Part I of this 
book. In essence, the main message and findings of this chapter are as follows. 
For a monocomponent linear FM signal, the WVD is optimal for energy concentra- 
tion about the IF and for unbiased estimation of the IF. If a signal has nonlinear 
frequency modulation and/or multiple components, the WVD suffers from inner ar- 
tifacts and/or outer artifacts (cross-terms), respectively; in either case, some form 
of reduced interference quadratic TFD (RID) is to be preferred over the WVD. The 
design of RIDs is best undertaken by designing the desired kernel filter in the am- 
biguity domain, and using Fourier transforms to see the effects in the time-lag and 
time-frequency domains. To be a useful tool for practical applications, quadratic 
TFDs are expected to be real, to satisfy the global and local energy requirements, 
and to resolve signal components while reflecting the components’ IF laws through 
the peaks of their dominant ridges in the (t, /) plane. Several RIDs were designed 
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Table 3.3.3: Properties of the TFDs whose kernels are listed in Table 3.3.2. The window w(t) is 
assumed to be real and even. An exclamation (!) means that the property is always satisfied, while an 
asterisk (*) means that the property is satisfied subject to normalization of the window (for tm) or the 
value of the parameter (for TS and fs). Comments on the kernel type are in parentheses (thus). 



Distribution 


PROPERTY 

RE TM FM IF TD TS FS 


Wigner-Ville 


! I I I ! I I 


Levin (product kernel) 


!!!!!!! 


Born-Jordan (product kernel) 


! ! ! ! ! * * 


Modified B (LI kernel) 


! ! ! ! 


w-WVD (DI kernel) 


! * ! ! 


w- Levin 


! * ! ! 


ZAM 


| * | * 


Rihaczek (product kernel) 


! ! ! ! 


w-Rihaczek 


* i 


Page 


! ! ! ! ! 


Choi- Williams (product kernel) 


! ! ! ! ! 


B (separable kernel) 


I 


Spectrogram 


I 



using simple separable kernels, demonstrating the procedure to construct quadratic 
TFDs that meet the above requirements. Special-purpose quadratic TFDs can be 
easily designed to meet the specifications of particular applications (as we shall see 
in several articles in the following chapters). In general, however, the WVD and 
various RIDs are the most useful TFDs; the spectrogram, which has been widely 
used, is at best subsumed by quadratic RIDs (of which it is a special case), and at 
worst made obsolete by them. 

The remaining four Parts of the book elaborate on these issues, discuss advanced 
design methods for TFDs, and present a wide selection of methodologies, algorithms, 
and applications that demonstrate how to use time-frequency methods in practice. 
In particular, questions and issues such as how to select a TFD for a particular ap- 
plication, how to implement it and how best to apply it are covered. These chapters 
constituting the remaining four Parts include a number of Articles which tend to 
cover more advanced and detailed material, complementing and supplementing the 
tutorial introduction of Chapter 1. 

Part II of the book gives more details on some fundamental topics of TFSAP 
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such as TFD design and signal analysis in the (t, /) plane. It includes two chapters 
(4 and 5). Chapter 4 presents some advanced concepts for TF signal analysis, TF 
signal processing and TF system analysis. Chapter 5 presents a number of methods 
for designing TFDs, including a treatment of the ambiguity function. 

Part III of the book describes specialized techniques used in implementation, 
measurement and enhancement of TFDs. It includes three chapters (6, 7 and 8). 
Chapter 6 deals with the implementation and realization of TFDs; in particular, the 
formulation of discrete-time quadratic TFDs is presented and computation issues are 
described. Chapter 7 presents quality measures for TFDs and methods for perfor- 
mance enhancement. Chapter 8 describes methods and algorithms for multi-sensor 
and time-space processing used in applications such as sonar and telecommunica- 
tions. 

Part IV presents the key statistical techniques for TFSAP of random signals. It 
includes four chapters (9 to 12). Chapter 9 presents time- frequency methods for ran- 
dom processes and noise analysis; Chapter 10 describes methods for instantaneous 
frequency estimation; Chapter 11 deals with the field of time- frequency synthesis 
and filtering, including time-varying filter design; and Chapter 12 presents time- 
frequency methods for signal detection, classification and estimation. 

Part V describes a representative selection of TFSAP applications encompass- 
ing a wide range of fields and industries. It includes three chapters (13 to 15). 
Chapter 13 presents time-frequency applications in telecommunications; Chapter 14 
describes time-frequency methods in radar, sonar and acoustics; and chapter 15 de- 
tails a number of time-frequency methods for diagnosis and monitoring used in a 
wide range of diverse applications. 

A detailed description of the contents of each chapter can be found in the table 
of contents of the book and in the introduction to each chapter. For a more detailed 
search of a topic needed, e.g. for the explanation of an advanced concept, the reader 
is referred to the detailed index at the end of the book. 

For various reasons such as space limitations and publication delays, a num- 
ber of recent references could not be included in the relevant chapters. They are 
sufficiently important to be listed here for the sake of completeness. They include 
new theoretical developments on a range of topics discussed elsewhere in this book 
such as TFDs with complex argument [25], rotated (t, f) kernel filters [26], Gabor 
analysis [27], detection [28], spectrogram segmentation [29], time- frequency plane 
decomposition in atoms [30], the issue of positivity [31,32], IF estimation [33,34], 
feature extraction [35] , an illustration of the concept of BT product [36] , polyspec- 
tra [37] and Volterra series [38]. 
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